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In this work we present a detailed computation of the spectra of UHECR in the top-down scenario. We compare 
the spectra of hadrons obtained by two different methods in QCD and supersymmetric (SUSY) QCD with large 
primary energies yfs up to 10 16 GeV. The two methods discussed are a Monte Carlo (MC) simulation and the 
evolution of the hadron fragmentation functions as described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equations. The hadron spectra obtained by the two methods agree fairly well in the interesting energy 
range 10 _5 Mx < E < O.SMx (Mx is the energy scale of the process Mx > 10 12 GeV). We have also computed 
the spectra of photons, neutrinos and nucleons obtaining a good agreement with other published results. The 
consistency of the spectra computed by different methods allows us to consider the spectral shape as a signature 
of the production model for UHECR, such as the decay of super heavy relic particles or topological defects. 



1. Introduction 

Ultra High Energy Cosmic Rays (UHECRs) are 
still an open problem in astro-particle physics. 
The 11 Akeno Grand Air Shower Array (AG AS A) 
events with energy larger than 10 20 eV pQ con- 
tradict the expected suppression of the UHECR 
spectrum due to the interaction with the Cos- 
mic Microwave Background (CMB) radiation, the 
Graisen-Zatsepin-Kuzmin (GZK) cut-off On 
the other hand the HiRes data seems to be con- 
sistent with the GZK cut-off picture pQ. If the 
UHECR primaries are protons and if they propa- 
gate rectilinearly, as the claimed correlation with 
BL-Lacs at energy 4 — 8 x 10 19 eV implies, than 
their sources must be seen in the direction of 
the highest energies events with energies up to 
2 - 3 x 10 20 eV detected by HiRes, Fly's Eye and 
AGASA pp. At these energies the proton atten- 
uation length is only about 20 — 30 Mpc and no 
counterparts in any frequency band was observed 
in the direction of these UHECR events. This is a 
strong indication that CR particles with energies 
larger than 10 20 eV may have a different origin 
from those with lower energies. 

Another important point of this discussion is 
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related to the low energy part of the CR spec- 
trum, there are, infact, strong evidences that 
these CR with energies 1 x 10 18 eV < E > 7 - 
8 x 10 19 eV are extragalactic protons, most prob- 
ably from Active Galactic Nuclei (AGNs). This 
statement is based on some robust experimen- 
tal and theoretical evidences: (i) extensive air 
shower (EAS) data confirm protons as primaries, 
(ii) the dip seen with great accuracy in the data 
of AGASA, HiRes, Fly's Eye and Yakutsk, is a 
strong signature of the propagation of UHE pro- 
tons in the extragalactic space, (iii) the begin- 
ning of the GZK cut-off seen in the spectra of 
AGASA and HiRes. Excluding the correlation 
of UHECR with BL Lacs from the analysis be- 
comes also possible the propagation of protons 
in very strong magnetic fields. Nevertheless, also 
in this case, the lack of a nearby source in the 
direction of the highest energy events (e.g. at 
£~ 3 x 10 20 eV) remains a problem. In fact, for 
very strong field strengths B ~ 1 nG: the deflec- 
tion angle, 8 ~ l at t/r H = 3.7°S nG given by the 
attenuation length Z at t and the Larmor radius r# , 
is small and sources should be seen. 

Many ideas have been put forward aiming to 
explain the observed superGZK (E > (6 — 8) x 
10 19 eV) events: strongly interacting neutrinos 
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and new light hadrons as unabsorbed signal carri- 
ers, ^-bursts, Lorentz-invariance violation, Topo- 
logical Defects (TD) and Superheavy Dark Mat- 
ter (SHDM) (see [3j for reviews). The two last 
models listed above, that represent the most 
promising top-down models, share a common fea- 
ture: UHE particles are produced in the decay 
of superheavy (SH) particles or in their annihila- 
tion. In the case of TD they are unstable and in 
the case of SHDM long-lived particles. We shall 
call them collectively X particles. Annihilation 
takes place in the case of monopolonia, necklaces 
4 and SHDM particles within some special mod- 
els. From the point of view of elementary parti- 
cle physics all these processes proceed in a way 
similar to e + e~ annihilation into hadrons: two 
or more off-mass-shell quarks and gluons are pro- 
duced and they initiate QCD cascades. Finally 
the partons are hadronized at the confinement ra- 
dius. Most of the hadrons in the final state are 
pions and thus the typical prediction of all these 
models is the dominance of photons at the high- 
est energies E > (6 — 8) x 10 19 eV. Let us now 
concentrate our attention to the computation of 
the UHECR spectrum produced in the decay of 
X particles. 

The spectrum of hadrons produced in the de- 
cay/annihilation of X particles is another signa- 
ture of models with superheavy X particles. The 
mass of the decaying particle, Mx, or the en- 
ergy of annihilation \/s, is in the range 10 13 - 
10 16 GeV. The existing QCD MC codes become 
numerically unstable at much smaller energies, 
e.g., at Mx ~ 10 7 GeV. Moreover, the comput- 
ing time increases rapidly going to larger energies. 
In this work we will review our results obtained, 
in the computation of the top-down spectrum of 
UHECR, using two different computational tech- 
niques: one based on a MC computation scheme 
[5] and the other based on the DGLAP evolu- 
tion equations [5]. In both cases SUSY is in- 
cluded in the computation. Monte Carlo simu- 
lations are the most physical approach for high 
energy calculations which allow to incorporate 
many important physical features as the pres- 
ence of SUSY partons in the cascade and coher- 
ent branching. The perturbative part of our MC 
simulation scheme is similar to other existing MC 



codes it also includes in a standard way SUSY 
and hence is reliable. For the non-perturbative 
hadronization part an original phenomenological 
approach is used in Ref. 5 . The fragmenta- 
tion of a parton i into an hadron h is expressed 
through perturbative fragmentation function of 
partons D\ (x, Mx), that represents the probabil- 
ity of fragmentation of a parton i into a parton 
j with momentum fraction x — 2p/Mx, convo- 
luted with the hadronization functions fj(x, Qo) 
at scale Qo, that is understood as the fragmen- 
tation function of the parton i into the hadron 
h at the hadronization scale Qo ~ 1.4 GeV 
To obtain the fragmentation functions of hadrons 
one has: 

D?(x,M x )=J2 / -Dl(-,Mx)f?(z,Qo)(l) 

where the hadronization functions do not de- 
pend on the scale Mx- This important property 
of hadronization functions allows us to calculate 
fi(x, Qo) from available LEP data, D'l(x, Mx) at 
the scale Mx — Mz, and then to use it for the cal- 
culation of fragmentation functions D^(x, Mx) at 
any arbitrary scale Mx- Our approach reduces 
the computing time compared to usual MC sim- 
ulations and allows a fast calculation of hadron 
spectra for large Mx up to Mgut ■ 

The perturbative part of the MC simulation in 
Rcf. 5 includes standard features such as angu- 
lar ordering, which provides the coherent branch- 
ing and the correct Sudakov form factors, as well 
as SUSY partons. Taking into account SUSY 
partons results only in small corrections to the 
production of hadrons, and therefore a simplified 
spectrum of SUSY masses works with good ac- 
curacy. The weak influence of supersymmetry is 
explained by the decay of SUSY partons, when 
the scale of the perturbative cascade reaches the 
SUSY scale Ql USY ~ 1 TeV 2 . Most of the en- 
ergy of SUSY partons remains in the cascade in 
the form of energy of ordinary partons, left af- 
ter the decay of SUSY partons. The qualitatively 
new effect caused by supersymmetry is the effec- 
tive production of the Lightest Supersymmetric 
Particles (LSP), which could be neutralinos or 
gluinos. The fragmentation functions D^(x, Mx) 
at a high scale Mx can be calculated also evolv- 
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ing them from a low scale, e.g. M x = Mz, where 
they are known experimentally or with great ac- 
curacy using the MC scheme. This evolution 
is described by the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equation 6 which can 
be written as 

dtD h = J- ^P tJ (z) ® D*(x/z,t) , (2) 

*■ — ' ITT 

j 

where t = ln(s/so) is the scale, ® denotes the 
convolution / ® g = J z dx/xf(x)g(x/z), and 
is the splitting function which describes the emis- 
sion of parton j by parton i. Apart from the ex- 
perimentally rather well determined quark frag- 
mentation function D^(x, Mz), also the gluon 
fragmentation function Dg(x, Mz) is needed for 
the evolution of Eq. (J2J). The gluon FF can be 
taken either from MC simulations or from fits to 
experimental data, in particular to the longitudi- 
nal polarized e + e~ annihilation cross-section and 
three-jet events. 

The first application of the DGLAP method for 
the calculation of hadron spectra from decaying 
superheavy particles has been made in Refs. [7J. 
The most detailed calculations have been per- 
formed by Barbot and Drees [7] , where more than 
30 different particles were allowed to be cascad- 
ing and the mass spectrum of the SUSY parti- 
cles was taken into account. Although at Mz, 
which is normally the initial scale in the DGLAP 
method, the fragmentation functions for super- 
symmetric partons are identically zero, they can 
be calculated at larger scales t: SUSY partons 
are produced above their mass threshold, when 
their splitting functions are included in Eq. (j2J). 
In [5] we proved that this method is correct. Also, 
the LSP spectrum can be computed within the 
DGLAP approach [7J. 

In this paper we shall study the agreement of 
the two methods: MC and DGLAP equations 
for the calculation of spectra produced in the de- 
cay or annihilation of superheavy particles. We 
shall also compare the results obtained by differ- 
ent groups comparing the calculated spectra with 
recent ones measured by UHECR experiments. 

The paper is organized as follows: in Section 
2 we will compare the two computation schemes 



described, referring the reader to our paper for 
more details. Photon, neutrino and proton spec- 
tra, needed for UHECR calculations, are com- 
puted in Section 3 and compared with the spectra 
obtained in [?J. In Section 4 we will discus the 
consequences of our results for models of SHDM 
and TD in the explanation of the UHECR spec- 
trum. We will conclude in section 5. 

2. MC and DGLAP comparison 

In this Section we shall compare the hadron 
spectra computed by the two methods discussed 
above, MC and DGLAP. In Figure □ we plot 
the FFs D£(x, M x ) calculated by SUSY MC and 
SUSY DGLAP methods for M x = 1 X 10 16 GeV 
and i = q. In the DGLAP method the SUSY FFs 
have been evolved from the ones obtained with 
the SUSY MC at the scale V? = 10M S usy ~ 
10 TeV. One can see the good agreement be- 
tween DGLAP (solid curve) and MC (dotted 
curve). This good agreement holds also for other 
(lower) scales Mx and for other initial partons, 
e.g. gluon, squark or gluino. One can see that the 
MC and DGLAP spectra slightly differ at very 
low x and have a more pronounced disagreement 
at large values of x. The discrepancy at low x 
is due to coherent branching, that is included in 
the MC scheme while it cannot be embedded in 
the DGLAP one. At large x, the calculations by 
both methods suffer from uncertainties, particu- 
larly the MC simulation. In this region the results 
are sensitive to the details of the hadronization 
scheme (see, e.g., the problem of HERWIG 
with the overproduction of protons at large x). 
When one does not have the initial SUSY FFs 
from a MC simulation, the question arises how 
to proceed. As was first suggested by Rubin [7], 
the initial FFs can be taken as the ones for ordi- 
nary QCD at the low scale y/s — Mz, while the 
production of SUSY partons is included in the 
splitting functions assuming threshold behavior 
at Msusy ~ 1 TeV. We can check this assump- 
tion computing the SUSY FF in both ways. In 
Figure □ we present the SUSY FFs D%(x,M x ) 
for i = q and M x = 1 x 10 16 GeV, evolved from 
the initial scale yfs = Mz (dashed curve). The 
good agreement between the two DGLAP curves 
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Figure 1. Comparison of SUSY DGLAP and 
SUSY MC fragmentation functions for Mx — 
1 x 10 16 GeV with quark as a primary parton. 
SUSY DGLAP FFs are calculated for 10 TeV as 
the starting scale (solid line) and for Mz (broken 
line). SUSY MC FF is shown by dotted line. 



Figure 2. Comparison of nucleon spectra from 
present work computed with DGLAP equation 
(solid line), from Barbot & Drees |Z| (dashed line) 
and Sarkar & Toldra [7] (dotted line). All three 
calculations are performed with quark as initial 
parton for M x = 1 X 10 16 GeV. 



proves the validity of the assumption made. 

3. Photon, neutrino and nucleon spectra 

The spectra of photons, neutrinos and nucleons 
produced by the decay of superheavy particles are 
of practical interest in high energy astrophysics. 
These spectra Df(x, Mx) with a = 7, v, N can be 
also considered as FFs. Because the dependence 
on the type i of the primary parton is weak, we 
shall omit the index i from now on, keeping a as 
subscript. 

Till now we concentrated our discussion on the 
total number of hadrons (a = h) described by the 
FF Dh(x,Mx), but in fact we have performed 
similar calculations separately for charged pions 
and protons+antiprotons. The procedure of the 
calculations is identical to that already described 
for the DGLAP and MC computation schemes. 
For charged pions and protons+antiprotons we 
used experimental data from Refs. [Sj. Below we 
shall present results of our SUSY MC simulations 
in terms of FFs for all pions D n , all nucleons 
Dn and all hadrons £>/,.. We introduce the ra- 
tios £n(x) and e ff (x) as: Djq(x) — ejsr(x)Dh(x) 



and D n (x) = e 7T (x)Dh(x). The spectra of pions 
and nucleons at large Mx have approximately 
the same shape as the hadron spectra, and one 
can use in this case e T = 0.73 ± 0.03 and £jv = 
0.12 ± 0.02 taking into account the errors in 
the experimental data We can calculate now 
the spectra of photons and neutrinos produced 
by the decays of pions neglecting the small con- 
tribution (0.15 ± 0.04) of K, D, A and other par- 
ticles. Including these particles affects stronger 
neutrinos than nucleons and photons, which are 
the main topic of this Section. 

The normalized photon and neutrino spectrum 
from the decay of one X particle at rest can be 
computed using the pion and nucleon FF fol- 
lowing the recipe given in (Jy. We shall com- 
pare our photon spectra with those calculated by 
the DGLAP method in |7j. The photon spec- 
trum is most interesting to compare, because it is 
straightforwardly related to the hadron spectrum 
which is the basic physical quantity. Moreover, 
the photon spectrum is the dominant component 
of radiation produced by superheavy particles. 

To be precise, we compare the FF 
Dl' N ' v {x,M x ) at M x = lx 10 16 GeV. Figure H 
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Figure 3. Comparison of SHDM prediction with 
the AGASA data. The calculated spectrum of 
SHDM photons is shown by dotted curves for 
two different normalizations. The dashed curve 
gives the spectrum of extragalactic protons from 
uniformly distributed astrophysical sources. The 
sum of these two spectra is shown by the thick 
curves. The x 2 values are given of the compar- 
ison of these curves with experimental data for 
E>4x 10 19 eV. 



that is referred to nucleons, demonstrates good 
agreement between our spectrum and those from 
at x < 0.3. As it was mentioned above, the 
disagreement at large x is not surprising. Apart 
from Dq(x, Mz) taken directly from the exper- 
iments, both calculations use the much more 
uncertain Dg(x,Q 2 ). In our case, D^{x,Q 2 ) is 
taken from our MC simulation jS], in the case of 
Barbot & Drees in [7] from the fit performed in 
Ref. 9 . In both cases, rather large uncertainties 
exist at large x jOJ. The agreement between the 
three curves as presented in Figure |21 is good. 
The same conclusion holds for the comparison of 
photons and neutrino spectra [S]- 

4. UHECR from SuperHeavy particles and 
Topological Defects 

As follows from the previous section, the accu- 
racy of spectrum calculations has reached such a 
level that one can consider the spectral shape as a 



Figure 4. Diffuse spectra from necklaces. The 
upper curve shows neutrino flux, the middle - pro- 
ton flux, and two lower curves - photon fluxes 
for two cases of absorption. The thick continu- 
ous curve gives the sum of the proton and higher 
photon flux. 



signature of the model. The predicted spectrum 
is approximately cx dE/E 19 in the region of x at 
interest. Another interesting feature of these new 
calculations is a decrease of the ratio of photons 
to nucleons, 7/-/V, in the generation spectrum. At 
i~lx 10~ 3 this ratio is characterized by a value 
of 2 - 3 only The decrease of the j/N ratio 
is caused by a decrease of the number of pions 
in the new calculations. This result has an im- 
portant impact for SHDM and topological defect 
models because the fraction of nucleons in the pri- 
mary radiation increases. However, in both mod- 
els photons dominate (i.e. their fraction becomes 
> 50%) at E > (7 - 8) x 10 19 cV. In this Section 
we shall consider: UHECR from superheavy dark 
matter (SHDM) [HI and topological defects (TD) 
[TT] . Production of SHDM particles naturally oc- 
curs in the time- varying gravitational field of the 
expanding universe at the post-inflationary stage. 
The relic density of these particles is mainly de- 
termined (at fixed reheating temperature and in- 
flaton mass) by their mass Mx- The range of 
practical interest is (3 — 10) x 10 13 GeV, at larger 
masses the SHDM is a subdominant component 
of the DM. SHDM is accumulated in the Galac- 
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tic halo with the overdensity S — = n PpM — , 

•> Px siCDMPcr ' 

where w 0.3 GeV/cm 3 , /9 cr = 1.88 x 10" 29 /i 2 
g/cm 3 and Ocdm^ 2 = 0.135 Q2|. With these 
numbers, S w 2.1 x 10 5 . Because of this large 
overdensity, UHECRs from SHDM have no GZK 
cutoff. 

dumpiness of SHDM in the halo can provide 
the observed small-angle clustering. The ratio 
rx = Six (to /t~x) of relic abundance fix and life- 
time tx of the X particle is fixed by the observed 
UHECR flux as rx ~ 10~ n . The numerical value 
of rx is theoretically calculable as soon as a spe- 
cific particle physics and cosmological model is 
fixed. In the most interesting case of gravitational 
production of X particles, their present abun- 
dance is determined by their mass Mx and the 
reheating temperature Tr. Choosing a specific 
particle physics model one can fix also the life- 
time of the X particle. There exist many models 
in which SH particles can be quasi-stable with 
lifetime tx ^S> 10 10 yr. The measurement of the 
UHECR flux, and thereby of rx , selects from the 
three-dimensional parameter space (Mx,Tr,tx) 
a two-dimensional subspace compatible with the 
SHDM hypothesis. 

In Figure |3| we have performed a fit to the 
AGASA data using the photon flux from the 
SHDM model and the proton flux from uniformly 
distributed astrophysical sources. For the latter 
we have used the non-evolutionary model. The 
photon flux is normalized to provide the best fit 
to the AGASA data at E > 4 x 10 19 eV. The 
fits are shown in Figure |21 with x 2 /d-°-f- indi- 
cated there. One can see from the fits in Fig- 
ure [3J that the SHDM model with the new spec- 
tra can explain only the excess of AGASA events 
at E > 1 x 10 20 eV: depending on the SHDM 
spectrum normalization and the details of the cal- 
culations for the extragalactic protons, the flux 
from SHDM decays becomes dominant only above 
(6-8) x 10 19 eV. 

Topological Defects (for review see [3] and ref- 
erence therein) can naturally produce UHE par- 
ticles. The following TD have been discussed 
as potential sources of UHE particles: supercon- 
ducting strings, ordinary strings, monopolonium 
(bound monopole-antimonopole pair), monopolo- 



nia (monopole-antimonopole pairs connected by 
a string), networks of monopoles connected by 
strings, vortons and necklaces (see Ref. for a 
review and references). Monopolonia and vortons 
are clustering in the Galactic halo and their ob- 
servational signatures for UHECR are identical 
to SHDM. However the friction of monopolonia 
in cosmic plasma results in monopolonium life- 
time much shorter than the age of the universe. 
Of all other TD which are not clustering in the 
Galactic halo, the most favorable for UHECR are 
necklaces. Their main phenomenological advan- 
tage is a small separation which ensures the ar- 
rival of highest energy particles to our Galaxy. 
We shall calculate here the flux of UHECR from 
necklaces. 

Necklaces are hybrid TD produced in the sym- 
metry breaking pattern G — > H xU(l) — > H x 
Z2- At the first symmetry breaking monopoles 
are produced, at the second one each (anti-) 
monopole get attached to two strings. This 
system resembles ordinary cosmic strings with 
monopoles playing the role of beads. Necklaces 
exist as the long strings and loops. The sym- 
metry breaking scales of the two phase transi- 
tions, rjm and n s , are the main parameters of the 
necklaces. They determine the monopole mass, 
m ~ Anr/m/e, and the mass of the string per 
unit length /i ~ ^Vs- The evolution of neck- 
laces is governed by the ratio r ~ m/ixd, where 
d is the average separation of a monopole and 
antimonopolc along the string. As it is argued 
in Ref. 0], necklaces evolve towards configura- 
tion with r> 1. Monopoles and antimonopolcs 
trapped in the necklaces inevitably annihilate in 
the end, producing heavy Higgs and gauge bosons 
(X particles) and then hadrons. The rate of X 
particles production in the universe can be esti- 

mated as @] fix ~ t^M x ' wnere * ^ s the cosmo- 
logical time. 

The photons and electrons from pion de- 
cays initiate e-m cascades and the cascade 
energy density can be calculated as uj caa = 

y^ 2 ^Io°wTnw = 2^ r2 #> where z is the 

redshift and f n ~ 1 is the fraction of the total 
energy release transferred to the cascade. The 
parameters of the necklace model for UHECR 
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are restricted by the EGRET observations [T3*| 
of the diffuse gamma-ray flux. This flux is pro- 
duced by UHE electrons and photons from neck- 
laces due to e-m cascades initiated in collisions 
with CMB photons. In the range of the EGRET 
observations, 10 2 — 10 5 MeV, the predicted spec- 
trum is cx E~ a with a = 2 [Tij. The EGRET 
observations determined the spectral index as 
a = 2.10 ± 0.03 and the energy density of ra- 
diation as w bs ~ 4 x 10~ 6 eV/cm 3 . The cascade 
limit consists in the bound w cas < L0 o bs- 

According to the recent calculations, the Galac- 
tic contribution of gamma rays to the EGRET 
observations is larger than estimated earlier, and 
the extragalactic gamma-ray spectrum is not de- 
scribed by a power-law with a = 2.1. In this 
case, the limit on the cascade radiation with 
a — 2 is more restrictive and is given by w cas < 
2 x 10 _6 eV/cm 3 ; we shall use this limit in fur- 
ther estimates. Using Lu cas with f n = 1 and 
to = 13.7 Gyr ^2] we obtain from the limit on 
the cascade radiation r 2 \i < 8.9 x 10 27 GeV 2 . 

The important and unique feature of this TD 
is the small separation D between necklaces. It 
is given by D ~ r -1 / 2 io 4 . Since r 2 /i is lim- 
ited by e-m cascade radiation we can obtain a 
lower limit on the separation between necklaces 

as D ~ (iffc) V4 *o > 10 WW 6 GeV 2 ) 1 / 4 kpc, 
this small distance is a unique property of neck- 
laces allowing the unabsorbed arrival of parti- 
cles with the highest energies. The fluxes of 
UHECR from necklaces are shown in Figure 
We used in the calculations r 2 [i = 4.7 x 10 27 GeV 2 
which corresponds to uj cas = 1.1 x 10 -6 eV/cm 3 , 
i.e. twice less than allowed by the bound on 
CL> cas . The mass of the X particles produced by 
monopole-antimonopole annihilations is taken as 
Mx = 1 x 10 14 GeV. ,/From Figure 0] one can 
see that the necklace model for UHECR can ex- 
plain only the highest energy part of the spec- 
trum, with the AGASA excess somewhat above 
the prediction. This is the direct consequence of 
the new spectrum of particles in X decays ob- 
tained in this work. Thus UHE particles from 
necklaces can serve only as an additional compo- 
nent in the observed UHECR flux. 



5. Conclusions 

In this paper we have compared the MC and 
DGLAP methods for the calculation of hadron 
spectra produced by the decay (or annihilation) 
of superheavy X particles with masses up to 
Mqut ~ 1 x 10 16 GeV. We found an excellent 
agreement of these two methods. The calcula- 
tions have been performed both for ordinary QCD 
and SUSY QCD. The inclusion of SUSY partons 
in the development of the cascade results only in 
small corrections, and it justifies our computation 
scheme with a single mass scale Msusy [E]< 

In comparison to the DGLAP method, the MC 
simulation has the advantage of including coher- 
ent branching. It allows reliable calculations at 
very small x. The Gaussian peak, the signature 
of the QCD spectrum, cannot be obtained using 
the DGLAP equations. We have calculated the 
all-hadron spectra, as well as spectra of charged 
pions and nucleons, using the DGLAP equations. 

Our nucleon spectrum agrees well with that of 
Refs. 0. We compared also our spectrum of pho- 
tons with the calculations of Ref. [Jj . The com- 
parison of the photon spectra is interesting, be- 
cause of physical reasons (photons can be observ- 
able particles), and because the photon spectra 
are connected directly with the hadron spectra. 
The spectra are in good agreement. We conclude 
that all calculations are in a good agreement es- 
pecially in the most interesting low x regime and 
the predicted shape of the generation spectrum 
(cx dE/E 1 - 9 ) can be considered as a signature of 
models with decaying (annihilating) superheavy 
particles. 

The predicted spectrum of SHDM model can- 
not fit the observed UHECR spectrum at 1 x 
10 18 eV < E < (6 - 8) x 10 19 eV. Only events 
at E > (6 - 8) x 10 19 eV, and most notably the 
AGASA excess at these energies, can be explained 
in this model. The robust prediction of this model 
is photon dominance. In present calculations this 
excess diminishes to j/N ~ 2—3 According to 
the recent calculations the muon content of pho- 
ton induced EAS at E > 1 x 10 20 eV is high, but 
lower by a factor 5-10 than in hadronic showers. 
The muon content of EAS at E > 1 x 10 20 eV has 
been recently measured in AGASA pQ . The mea- 
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sured value is the muon density at the distance 
1000 m from the shower core, p p (1000). From 
11 events at E > 1 x 10 20 eV the muon density 
was measured in 6. In two of them with ener- 
gies about 1 x 10 20 eV, is almost twice higher 
than predicted for gamma- induced EAS. Taking 
into account the contribution of extragalactic pro- 
tons at this energy, the ratio j/p predicted by 
the SHDM model is 1.2 - 1.4. It is lower than 
the upper limit 7/p < 2 obtained by AG AS A at 
E = 3 x 10 19 eV on the basis of a much larger 
statistics. The muon content of the remaining 
4 EAS marginally agrees with that predicted for 
gamma-induced showers. The contribution of ex- 
tragalactic protons for these events is negligible, 
and the fraction of protons in the total flux can be 
estimated as 0.25 < p/tot < 0.33. This fraction 
gives a considerable contribution to the probabil- 
ity of observing 4 showers with slightly increased 
muon content. Not excluding the SHDM model, 
the AG AS A events give no evidence in favor of it. 

The simultaneous observation of UHECR 
events in fluorescent light and with water 
Cherenkov detectors has a great potential to 
distinguish between photon and proton induced 
EAS. An anisotropy towards the direction of the 
Galactic Center is another signature of the SHDM 
model. Both kinds of informations from Auger £Q 
will be crucial for the SHDM model and other 
top-down scenarios. Topological defect models 
are another case when short-lived superheavy 
particle decays can produce UHECR. In Figure 
SI the spectra from necklaces are presented. One 
can see that at E > 1 X 10 20 eV photons dom- 
inate, and the discussion in the previous para- 
graph applies here too. In contrast to previous 
calculations, the agreement with observations is 
worse: necklaces can explain only the highest en- 
ergy part of the spectrum in Figure 01 with the 
AGASA excess somewhat above the prediction. 
In the other energy ranges, UHE particles from 
necklaces can provide only a subdominant com- 
ponent. Other TDs suffer even more problems. 
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